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0.  Summary 


We  will  find  the  probability  distribution  of  losses  from  reflections 
at  discontinuities  in  a  transmission  line,  on  the  assumption  that  the 
reflection  coefficients  at  the  discontinuities  have  random  phase,  and,  in 
addition,  random  magnitudes  associated  with  manufacturing  uncertainties 
(tolerances)* 

1.  Introduction 

Suppose  a  transmission  line  has  n  discontinuities  with  associated 

voltage  steuiding  wave  ratios  (VSWR’s)  <  r2  <  <  *, ,  <  A 

common  procedure  for  describing  the  cumulative  effect  of  reflections  at 

the  discontinuities  is  to  say  that  the  resulting  VSWR  will  be 

r-r^r_,,.r  when  the  mismatches  add  in  the  worst  phase,  or  r - 
J.4i^n  ^  n  X  d  n— X 

when  they  combine  in  the  best  phase,  (See,  for  example,  [4],  p.35f  6ll). 
This  gives  a  rather  large  intervad  which  is  of  limited  use  to  designers. 

Can  a  more  precise  description  of  the  overall  VSWR  be  given  in  probability 
terms?  In  [5],  Mullen  and  Pritchard  give  such  a  description,  based  on 
the  approximation  that  reflection  coefficients  are  additive.  We  will 
develop  a  procedure  based  on  combining  the  reflection  coefficients  in  a 
more  precise  way, 

2,  Combining  Reflection  Coefficients 

Let  Y21  reflection  coefficients  associated  with  n 

discontinuities  in  a  transmission  line.  We  aire  interested  in  the  overall 
reflection  coefficient  y.  assume  that  •••t  cm  not  be 
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given  exactly,  €uid  hence  that  each  ie  a  random  complex  number  vrith  a 
distribution  that  will  be  specified.  We  want  the  probability 
distribution  of  Yi  particular,  that  of  |yI  i  which  will  give 

an  indication  of  the  total  reflection  losses  in  the  system. 

Now  a  rough  approximation  to  y  is  Y  =  Yj^  +..•+  Y^^i  was 
considered  in  [5  3.  We  will  use  a  more  precise  relation  between  y 
and  the  find  the  probability  distribution  of  y  from  that  of 

the  relation  we  want  is 


This  is  equivalent  to  adding  the  shunt  admittances  of  the  discontinuities. 
See,  for*example,  Montgomery,  [2J,  p,71. 

If  we  let  these  generalized  admittances  be 

* 

n  =  y/(1  -  =  Y^AI  -  Y3^)|. ..»»!„  =  Y^O  -  Yq)» 

then 

n  =  +  ^2  +...+  \ 

and  Y  =  >l/(l  +  »l)»  will  assume  that  are  Independent 

random  conQ)lex  variables,  and  hence,  by  the  central  limit  theorem,  we 
can  expect  that  the  real  and  complex  p8u*ts  of  will  be  jointly 

normally  distributed.  We  then  must  find  the  distribution  of  jylt  vhen 
Y  =  VO  ’l)  distribution  of  n  is  the  bivariate  normal. 
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5 •  Distribi J lion^ of  the  Adjnittances 

Our  procedure  for  finding  the  distribution  of  the  overall  nti.wvxon 
coefficient  is  this:  we  are  given  n  reflection  coefficients 
^l’^2* * • • *^n’  assumed  to  be  independent  random  complex  numbers.  We  first 
convert  reflection  coefficients,  io  admittances,  according  to 

the  relation  -  Yj^).  We  then  add  the  admittances,  getting 

H  =  T]^.  By  the  central  limit  theorem,  r\  will  be  approx¬ 

imately  normally  distributed  in  the  complex  plane.  Then  we  convert  y\ 
back  to  an  overall  reflection  coefficient,  y,  according  to  the  relation 
y  =  Ji/(1  +  Tt). 
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we  may  take  x  and  y  to  be  independent  normal  random  variables  with 
means  and  variances  calculated  from  expression  (P)  in  the  appendix «  In 

particular  y  when  we  assume  that  each  is  uniformly  distributed  over 

(0,2n),  then  we  have 

• 

E(x)  =  E{y)  =  0 

4,  The  Distribution  of  lii/d  ^  t|)|  from  that  of  , 

We  have  y)  =  x  'f  iy,  where  x  and  y  are  independent  normal  random 
variables,  and  want  the  distribution  of  lylf  where  y  =  t/d  + 
Presumably,  |yI  takes  values  between  0  and  1,  euad  for  0  <  t  <  1, 

we  note  that  IyI  <  ^  only  if,  in  the  complex  plane,  the  distance 

from  y\  to  the  origin  is  less  than  t  times  the  distance  from  x\  to 
-  1  Oi,  In  terms  of  x  and  y,  this  condition  may  be  put  in  the  form 

pCIyI  <  t]  =  p[(x  -  tVd  -  t^))  +  <  (t/(i  - 

Thus  the  probability  that  IyI  will  be  less  than  t  la  the  bivariate 

normal  probability  measure  of  an  offset  circle.  The  standard  bivariate 

normal  probability  measure,  (a  =  o  =1),  of  a  circle  with  center 

X  y 

distauice  d  from  the  origin  and  radius  r,  designated  C(d,r), 

C(d,r)  =  PC(5  +  d)^  ♦  <  r^] 

is  a  well  known  function.  It  may  be  computed  by  the  formula 

C(d,r)  =  S  ^  F.  (r^) 
k=0  “ 
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where 


X  t.  d^/2,  and.  if  it  -  r^/2, 


a  1  -  e-'*  -  ne-'*  -  |j  a''* 


k 


Tables  of  C(d,r)  are  available  from  the  author,  [  1  ],  or  from  Rand^  [5], 


Our  results  may  be  summarized  as  follows: 

If  a  transmission  line  has  discontinuities  with  reflection  coefficients 
iu.  iu^  iu^ 

=  c^e  ,Y2  =  ^2®  »  where  the  c*8  and  u*8  are 

independent ,  each  with  some  distribution  on  the  interval  (0,1),  and 

each  u,  uniformly  distributed  on  the  interval  (0,2ti),  then  the  resulting 
reflection  coefficient >  Yt  related  to  by 


has  magnitude  with  probability  distribution 

(1)  PCIyI  <  t]  =  C[  - ]  . 

o(l  -  t^)  o(l  -  t"^) 

where 

0^  =  I  E  B[c^/tl  -  cj)] 

and  C(d,r)  is  the  standard  normal  probability  measure  of  a  circle  with 
radius  r,  center  distance  d  from  the  origin^ 

Graphs  of  the  distribution  function  and  density  function  of  IyIi  for 
0  =  .02,. 04,...,, 26,  are  drawn  in  figure  1. 
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We  have  seen  that  the  probability  distribution  of  the  magnitude^ 
of  the  resultant  reflection  coefficient,  on  the  assumption  that  the 
phases  of  the  component  reflection  coefficients  are  uniformly  distributed, 
is  given  by  (!)♦  It  depends  on  a  single  parameter,  o,  related  to  the 
individual  reflection  coefficients 

=  I  S  ECc^/tl  -  C^)]. 

It  is  therefore  important  that  we  determine  the  expected  value  of 
2  2 

•  c^/tl  -  Cj^)  for  each  k.  In  most  practical  problems,  c  will  be  given 
in  terms  of  a  tolerance,  e.g,,  c  =  .10  +  ,04.  How  can  we  determine 

E(c  /(I  -  c  )  when  all  we  know  is  that  c  is  given  in  a  form  such  as 

c  =  .10  +  ,04?  Fortunately,  we  can  get  good  upper  and  lower  bounds  on 
this  expectation  for  two  reasons s  The  function  x  /(l  -  x  )  is  convex, 
and  the  range  of  the  typical  c  is  small.  The  bounds  we  want  are  as 

follows:  If  c  is  given  by  c  =  c^  +  6,  where  c^  is  the  expected  value 

of  c,  then  no  matter  how  c  is  distributed  on  the  Interval 
Cq  -  6  <  c  <  Cq  -f  6,  we  have 


where  L  =  c^  -  6  is  the  lower,  and  U  =  c^  •••  6  the  upper,  tolerance 
value  on  c  .  In  the  example  above,  c  =  .10  ^  .04,  we  have 
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or 

2 

.0101  <  E( - 2__)  <  ^0118. 

1  -  c 

Thus  we  have  contained  our  unknown  expected  vailue  in  a  rather  small 
interval. 


6.  •  Examples 


Example  1.  *  Suppose  we  have  8  discontinuities ^  each  with  VSWH  rated 
1,1  +  ,04.  Conventional  theory  states  that  the  worst  mismatch  will  result 

Q 

in  a  VSWR  of  about  (1,1)  =  2,l4,  or  a  reflection  coefficient  with 

magnitude  (2,l4  -  1)/(2,14  +  1)  =  ,36,  Figure  2  gives  the  actual 

distribution  of  the  resultant  lylt  it  is  clear  that  only  very  rarely 

will  IyI  be  anywhere  near  ,36,  To  describe  the  possible  variation  of 

IyI  in  probability  terms,  we  first  convert  tolerances  on  each  VSWH  to 

tolerances  on  the  magnitude  of  the  corresponding  reflection  coefficient. 

If  r  is  in  the  range  1,06  <  r  <  I,l4  then  c  =  (1  -  r)/(l  +  r)  is  in 

the  range  ,029  <  c  <  .065,  which  we  write  as  c  =  ,048  +  ,0l8,  taking 

,048  BB  the  expected  value  of  c.  Then  we  assume  that  each  reflection 

coefficient  Yj^  =  ^  ^  randcwn  complex  number  with 

Independent,  vd.th  some  distribution  on  .048  ♦  ,0l8  and  uniform 

2  2 

on  (0,2n).  We  then  need  According  to  section  5,  we 

have  upper  and  lower  bounds: 


c^  ana  uj^ 


.048* 


1  -  .048 


.066^ 


1  -  .066* 


or 

.0023  <  ECc^/l  -  c^)]  <  .0026. 
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Thus  if  11  =  Ti^  +  ^2  +...+  Tij^  =  -  Yj^)  i  then  n  =  x  +  iy 

has  components  that  are  (approximately)  independent  normal  random 
variables  with  zero  means  and  common  variance 

^  =  |SE[C^/<1  -  c^)]  <|(8)(,0026)  =  .0104. 

We  take  c  -  ,1,  and  the  magnitude  of  y,  the  resulting  reflection 
coefficient ,  has  distribution 

p[|yI  <  t]  =  p[(x  -  t^/(i  -  t^))^  +  <  (t/(i  -  t^))^] 

or 

pcIyI  <  t]  =  cdot^/d  -  t^),iot/d  -  t^)). 

Figure  2  gives  the  distribution  and  density  function  of  lyl  for 
^  this  example.  It  is  clear  that  the  resulting  value  of  lyl  is  virtually 
certain  to  be  much  better  tham  the  pessimistic  value  of  ,36  obtained  by 
multiplying  the  VSWR’s  -  the  probability  that  M  <  .21  is  about  ,90, 
and  half  the  time  lyl  will  be  less  than  ,12, 


Figure  2.  Distribution  and  densitjj  function  of  |y|>  Example  !• 
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Example  2,  Suppose  we  have  15  discontinuities,  11  with 
c  -  .05  +  .02  and  4  with  c  =  ,10  +  .04.  Also,  suppose  an  additional, 
final,  discontinuity  at  the  antenna  has  complex  reflection  coefficient 
known  to  be  .5©*^^^^.  What  can  we  say  about  total  reflection  losses  in 
the  system?  For  c*s  in  the  range  ,05  +  ,02,  we  have 


<  * 


1  -  .05^  "  1  -  c 


1  -  .03  1  -  .Cf? 


.0025  <  E( - 2_)  <  ^0029  • 

1  -  c 

and  for  c*s  in  the  range  .10  +  .04,  we  found,  in  section  5| 

2 

.0101  <  E(— -— )  <  .0118. 


Thus  the  generalized  admittance,  + 

of  the  random  admittances  of  the  15  discontinuities,  has  the  form 
i|  =  X  +  iy,  where  x  and  y  are  independent,  normal,  zero  means,  common 
variance 

<  (11) (.5) (.0029)  +  4(.5)(.01l8)  =  .0396. 

We  thus  talce  0  =  V,0396  =  .2  . 


To  tf  =  X  +  iy  we  must  add  the  generalized  admittance  of  the  last 
discontinuity,  which  is  .5©*^^^^!  or,  in  rectangular 

coordinates,  -  ,191  >,651i. 

Thus  the  generalized  admittance  of  the  entire  I6  discontinuities  is 


T  =  X  -f  .191  +  (y  .651)i 


.  ^ 
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and  if  Y  =  V(1  +  '*’ )  is  the  overall  reflection  coefficient, 

PCIyI  <  t]  =  P[(x  +  .191  -  t^/(l  -  t^))^  +  (y  +  .651)^  <  (t/(l  - 
or 

<  t]  .  p[(i  .  ^ - .  (I  *  <  ( - L^)4]. 

®  ®  d(l  -  t^)  °  ®  ~  0(1  -  t^) 


The  last  expression  is  the  standard  normal  probability  measure  of  a  circle 

of  radius  t/a(l  -  t^),  center  distance  0*^[(.191 - ^ (*651)^]^^ 

•  •  1  -  t2 

from  the  origin,  i.e.^  with  o  =  .2, 

•  PCIyI  <  t]  =  C(d,r) 

with 

r  =  5t/(l  -  t^),  d  =  Me. 191  -  t^/(l  -  t^)]^  +  (.651)^}^^. 

Figure  3  gives  the  distribution  and  density  functions  of  IyI  for 
this  example. 


Figure  3.  Distribution  and  density  of  |y|»  Exanple  2. 
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APPENDIX 


Theorem  1. 

Let  Y  =  ce^^  be  a  complex  random  variable  in  polar  form^  where 
c  and  u  are  independent  random  variables,  and  c  has  some  dietributlon 
In  the  interval  (0,1). 

If  u  is  symmetrically  distributed  over  the  interval  (0,2n),  then  the 
real  and  imaginary  parts  of  tj  =  y/d  -  y)  are  uncorrelated. 

If  we  compute  the  real  and  imafrinary  parts  of  y\  directly,  we  have 

//,  V  .  c  cos  u  -  ic  sin  u 

=  y/(l  -  y)  =  X  +  iy  =  — ;; - ^  , 

l-2c  cos  u+c  l-2c  cos  u+c 


and  it  is  not  evident  that  x  and  y  are  uncorrelated,  fowever,  if  we 

2  3 

write  r\  =  y/d  -  y)  -  y  +  y  y^  then 

2  5  2  3 

(1)  II  =  (c  cos  u  +  c  cos  2u  +(T cos  3u  -f  i<c  sin  u  -f  c  sin  2u  +  c^sin  3u 

and  we  may  take  expectations  term  by  term  to  get 

2  3 

E(x)  =  E(c)E(co6  u)  +  E(c  )£(cos  2u)  +  £(c^)E(co8  3^)  +,,, 

(2)  E(y)  =  E(c)E(sin  u)  E(c^)E(sin  2u)  >  E(c^)E(sin  3u) 

E(xy)  =  22  E(c^'*‘°)E(co8  ru  sin  su), 
rs 

Now  if  u  is  symmetricailly  distributed  over  (0,2k)  then 

E(sin  ru)  =  E(sin  ru  cos  su)  =  0,  Hence  E(y)  =  0  and  £(xy)  =  0,  so 

that  X  and  y  are  uncorrelated. 
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Of  particular  interest  is  the  case  where  u  is  uniformly  distributed: 
Theorem  2. 

^  Y  =  ce^^,  where  c  and  u  are  independent,  c  has  some 
distributions  on  the  interval  (0,1),  and  u  is  uniformly  distributed 
on  the  interval  (0,2ti),  then  the  real  and  imaginary  parts  of 
=  y/(1  -  y)  =  X  +  iy  are  uncorrelated,  with  expected  values 

E(x)  =  E(y)  =  E(xy)  =  0 

and  variances 

2 

E(x^)  =  E(y^)  =  |eC— . 

The  proof  of  Theorem  1  applies,  except  that  we  must  compute  the 

variances  of  x  and  y.  Squaring  and  taking  expected  values  of  the  series 

expressions  for  x  or  y,  the  cross  product  terms  are  zero,  and  there 

2  2  1 

remains,  since  ECsin  ru)  =  E(cos  ru)  =  , 

E(x2)  =  EC/)  =  |[E(c^)  +  ECc'^)  +  E(o^)  +...] 
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Figure  1.  Distribution  and  density  functions  of  |yI»  ilie  magnitude  of 
the  overall  reflection  coefficient,  for  various  values  of  a. 


